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Abstract: The work is devoted to mathematical modeling of the drying process of potash 
fertilizers in fluidized bed apparatuses. The model is based on the material and thermal balance and 
takes into account the kinetics of the dehydration process of potassium chloride particles. The model is 
investigated using the tools of the Matlab mathematical computing package. The adeguacy of the 
proposed analytical model to real processes and phenomena has been verified. 


Keywords: Mathematical modeling and optimization, potash fertilizers, drying of potassium 
chloride in asuspended layer, virtual guality analyzer. 


There is a significant amount of work devoted to modeling the drying processes of bulk materials, 
in particular, drying in fluidized bed apparatuses. At the same time, numerous models describing the 
structure of the fluidized bed reflect only certain aspects of the phenomenon. Due to the complexity of 
the system and with a large number of internal and external parameters, these models are suitable for 
use as a basis for reflecting empirical correlations and dependencies [1-2]. The mathematical model is 
based on the material and energy balance of drying, using simplifications and empirical dependencies 
defined, for example, in [3]. The material and energy balance of dehydration is determined by the 
following eguations: 


M(1-U,)=(M-U)( 1-U,) (1) 
U=R(4); (2) 
qı 
qı = (Quen T 2, OI) +(m —1)C rt ; (3) 
sa a 
(UU 
TAR 
q=(Q+0,6t,, 1915 >) (5) 


F 
where U;.— final humidity, (mass fraction); U„— initial humidity, (mass fraction); U— the amount 
of evaporated moisture, (kg/s); g— the amount of heat, transmitted in a layer by a coolant per one kg of 
burned fuel, (J/k); qi— the amount of heat per evaporation of one kg of moisture, (J/kg); m— consumption 
of wet material corresponding to one kilogram of evaporated moisture, (ke/ke Barn); Cine — heat 
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capacity of the material (potassium chloride-KCl), (J/(kg:K); tea — fluidized bed temperature, (°C); f»— 


coolant temperature, (°C); Q- calorific value of gas, (Lorcixe), qucn — the heat of evaporation of water at 
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09C, (J/kg); n— coefficient of heat loss in the furnace and fluidized bed apparatus; M— the amount of 
wet material entering the drying, (kg/h); R— total fuel consumption, (kg/s). 

The gas temperature in the sublattice space , depending on the fuel gas consumption , is calculated 
by the empirical formula [4]: 


tr = 2920 -D / a+ 0,43. (6) 
In this case , the coefficient a is calculated by the formula 


a = L/(R-9,4) (7) 


where £ — coefficient that takes into account heat loss in the furnace; a — excess air ratio; L 
— hourly secondary air consumption, (nm?/s). 
As a result of the transformations we get 


(E 00,014(m-1)-C,,)-0,65 (2) 
Rm 


t 
baa = E 
09 
ty 
U Q 0.6, U i 
((—)(2, 01+ (m-1)-C,,)-0,6+ (DY -4 CN A rey A 

Rm lr tr RN 

F 1,2 
(—) 
tr 


The dependence of the final humidity on the temperature of the layer is expressed by the 
empirical formula: 
UA (9) 


where the coefficients A and k depend on a variety of process conditions, in particular on 


the ambient temperature, the granulometric composition of the product, and are determined 


1 
empirically and take the following values: A=0.02237, wt, ; k=-0,025, = 


C 


The model based on the material-thermal balance is implemented using the means of the 
MATLAB mathematical computing package. The operation of the model was tested on an array 
of historical data from the automated process control system archive obtained on functioning 
equipment. At the same time, sampling to determine the moisture content of the product in the 
laboratory was carried out at intervals of 600 seconds during the working shift. Figure 1 shows 
graphs of changes in the humidity of potassium chloride particles at the outlet of the fluidized 
bed dryer. In the graph shown in Fig.1, the series indicated by the dotted line is a sample of the 
values of the final moisture content of the Ux; product in the real process (Fig. 2, a). The series 
indicated by a solid line represents the result of the analytical model at the values of itsparameters 
indicated above, as well as the values of the crystallizate flow rate at the inlet and theinitial 


humidity indicated in Fig. 2, b. If there is a selective distribution law of random observations 
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(determined by the model), then the adeguacy of the analytical model can be assessed using the 
Pearson agreement criterion (y? - criteria) [4]. 


To apply the y7-criterion, the entire range of variation of a random variable in a sample of 
volume N is divided into r intervals. The number of intervals when calculating the Pearson 
consensus criterion is determined by the formula 


r =143,322-1g(N) (10) 
where N is the number of values in the sample, r is the number of intervals. 


The number of intervals r (for n — 45) is 6. Based on the maximum and minimum values of 
the residual humidity index in the sample under consideration and the specified number of 
intervals, the interval step of 0.02 wt is calculated. %. 


Change in humidity at the outlet 
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Fig.1. Change in the humidity of potassium chloride (KCI) particles at the outlet of the 
fluidized bed apparatus 
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Fig. 2. Results of the analytical model of the fluidized bed furnace 
Calculation of the freguencies of the results presented in Fig. (2,a), is given in Table.1 


Table.1 

Results of calculation of the analytical model 
Interval | Interval boundary, wt.% Values of the frequency of hits in 
number the interval 

lower upper settlement laboratory 

1 0,1 0,12 2 2 
2 0,12 0,14 8 10 
3 0,14 0,16 11 6 
4 0,16 0,18 11 9 
5 0,18 0,2 10 13 
6 0,2 0,22 3 5 

The calculated value of the y--criterion is calculated by the formula 

r 2 
pact yo (11) 


where ¥pac- calculated value of the criterion y7; O; - hit rate 

in the i-th interval of calculated values according to the model; E ;— the frequency of hitting the 
i-th interval of laboratory values; r is the number of intervals. 

The calculated value of the criterion y7 ras for the data in Table. 1 equals 6,5. According to the table 
of critical values of y ras , at the significance level a = 0.05 and the number of degrees of freedomë = r- 
3 - 6 — 3 — 3 , the value of y7 crete is equal to 7,8. At the same time, y ras < X? crete, i.e. the calculatedvalue 
of the criterion turned out to be less than the critical one.Thus, the hypothesis of eguality (agreement) 
of freguencies is not rejected and the model is recognized as adeguate. 

To assess the guality of the developed model, you can additionally use other indicators, for 
example: 

1. Correlation coefficient. Reflects the tightness of the linear relationship between the calculated 
values of the model and the experimental ones. 


p20 LED) (12) 
ND” GE yy 


2. Average approximation error 


yo» 


delta- ___>' -100 (13) 
N 


3. Standard approximation Error (RMSE 


Y T 4. Maximum approximation error 
yr I 

sko _ delta =s; 

= V N 
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max. delta = max(abc(y — 49 (15) 
Here: $4— calculated value according to the model; y;— experimental value; $ estimation of 


the mathematical expectation of the vector of calculated values; y — estimation of the mathematical 
expectation of the vector of experimental values; N— the number of values in the output vector of the 


model. 
The values of the quality indicators of the analytical model (12)-(15): 
correlation coefficient R........eeue..e.e ueu Y Y iN. 0,673 
average relative BITOF FF —— OY YO OOIRARŴH 0,006 
RMS mean square GETOL cnd eu Ry ddo 0,023 
PAK SL ON ETTO HCC EF FFAT TPFF 0,057 


Information for assessing the accuracy of the analytical model can also be provided by the errors 
given regarding the range of values of laboratory analyses in studies and the range of corresponding 
calculations based on the analytical model. The relative reduced error is calculated by the formula 


A 
Y — __ 10096 (16) 


n 


where Y — relative reduced error of the analytical model; A — average absolute measurement 
error; Xn — the normalizing value of the scale. 


The difference between the maximum and minimum values of the laboratory analysis 
vector, equal to 0.12 wt, % was used as X, and the difference between the maximum and 
minimum values of the vector of calculated values, equal to 0.1132 wt.%. 


The average absolute error of the analytical model is egual to: 


EP (17) 


The relative reduced error of the model for the range of values of laboratory analyses (16) 
is numerically equal to 14.590. The relative reduced error of the analytical model based on the 
range of calculations of the model (16) is numerically egual to 15.390. 
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